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Superfluidity, the ability of a fluid to move without dissipation, is one of the most spectacular 
manifestations of the quantum nature of matter. We explore here the possibility of superfluid motion 
of light. Controlling the speed of a light packet with respect to a defect, we demonstrate the presence 
of superfluidity and, above a critical velocity, its breakdown through the onset of a dissipative phase. 
We describe a possible experimental realization based on the transverse motion through an array of 
waveguides. These results open new perspectives in transport optimization. 
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Next year will bring the opportunity to celebrate the 
100th anniversary of the discovery of superconductiv- 
ity This remarkable property is often related to a 
more fundamental phenomenon, the Bose Einstein con- 
densation, where a single quantum state is occupied by 
a macroscopic number of particles. The bosons that 
condense may be coupled electrons that form Cooper 
pairs, as in superconducting metals 0, atoms [lb like 
in the original experiments in superfluid (SF) 4 He |4j, or 
molecules [5| . They may also be formed of more complex 
particles, like fermionic atom pairs @, or polaritons, a 
composite of a photon and an exciton Q. Superfluid- 
ity of polaritons in semiconductor cavities was explicitly 
tested recently 

There are different definitions of superfluidity, each one 
may emphasize a particuler physical aspect. Here super- 
fluidity means the existence of a finite critical velocity 
v c > below which the motion of the fluid is dissipation- 
less. A particularly simple way to experimentally imple- 
ment this test is by moving through the fluid an obsta- 
cle, or localized external potential. When the potential is 
weak, it has been shown long ago [Si that v c — c s , where 
c s is the speed of sound in the fluid [10( . Above the criti- 
cal velocity, superfluidity is broken and dissipative effects 
appear. 

It is important to note that a finite critical velocity is 
directly related to the presence of interactions between 
the bosons. The interactions control the long wavelength 
structure of the dispersion relation of the fluid. In a 
mean field approximation, weakly interacting bosons may 
be modeled, with good accuracy, by the Gross-Pitaevsky 
equation. The latter corresponds to a Schrodinger equa- 
tion with an additional nonlinear term that describes 
the interactions. In particular, this equation reproduces 
correctly the Bogoliubov dispersion relation mentioned 
above for the excitations of the interacting fluid. Interest- 
ingly, when the propagation of light is considered through 
a nonlinear medium of the Kerr-type uniform in one di- 
rection, in the paraxial approximation a similar equation 
is obtained for the slowly varying envelope of the optical 
field of a given wavenumber and frequency. The anal- 
ogy between the Gross-Pitaevsky equation and the light 



propagation in nonlinear media has been exploited in the 
past to test basic quantum effects with optics, like Bloch 
oscillations 11, 12 1 or Anderson localization [13j |. Be- 
cause of the similarities of the two equations, and since 
the Gross-Pitaevsky equation predicts SF motion, it is 
natural to push further the analogies and consider the 
possibility to observe a new state of light, e.g. superflu- 
idity in an optical nonlinear medium. Based on a self- 
defocusing refractive medium inside a Fabry-Perot cavity, 
an optical analog of a SF has indeed been proposed [14j . 
However, the results of the numerical simulations based 
on a transient regime were not conclusive. Moreover, no 
clear evidence of the existence of a SF critical velocity 
was provided. Therefore, the existence of photonic su- 
perfluidity in nonlinear media, as well as its experimental 
observation, remain open issues. 

Our purpose here is (i) to provide clear evidence of 
SF motion of light in a nonlinear medium as well as of 
its breakdown and (ii) propose an experiment that al- 
lows the observation of these effects. For simplicity, we 
will focus on the propagation of light in an effective one- 
dimensional array of waveguides. 

In these materials, light propagates in a medium where 
the refractive index has been spatially modulated. A 
typical set up consists of a periodic modulation of a two- 
dimensional layer, where an array of equally spaced iden- 
tical waveguides is formed (see Figfl]). Outside each of 
the waveguides the optical field intensity decreases expo- 
nentially. When the distance is such that the overlap be- 
tween the fields of neighbouring waveguides is small, the 
optical tunnelling between adjacent guides may simply 
be modelled by a hopping term. Light propagates along 
the guides in the longitudinal direction and hops from 
guide to guide in the transverse direction. Moreover, the 
width of each waveguide may be engineered in order to 
modify the energy of the local (quasi) bound state light 
mode, and Kerr materials may be used to include non- 
linear effects. Under these conditions, the optical field 
amplitude Ak of light at the k th lattice site (or waveg- 
uide) obeys the following discrete nonlinear Schrodinger 
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FIG. 1: (color online) Schematic view of an array of waveg- 
uides. The red arrows indicate the input beam and the output 
facet, where the light intensity distribution is measured. 



equation [12j (in the paraxial approximation) 
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where C is the tunneling rate between two adjacent sites, 
e k the on-site energy, and 7 > the strength of the 
self-dcfocusing nonlinearity of the medium. The left 
hand side describes the propagation along the longitu- 
dinal z-axis of the waveguide, and replaces time in the 
Schrodinger equation. We measure all distances in units 
of the incident light wavenumber, hence z, C, e k , and 
7|Afc| 2 are dimensionless. 

The possibility to engineer the different characteristics 
of the array make these photonic lattices unique in their 
ability to control the flow of light. A laser beam is shone 
on the input facet of an array of N waveguides, prop- 
agates across the photonic lattice, to finally reach the 
output facet where the intensity distribution is measured 
(as shown in FigH]). The length of the structure in the z- 
direction thus determines the time of propagation across 
the lattice. 

In order to test superfluidity of light, we are interested 
in the scattering properties of an incident pulse on a lo- 
calized defect. In the absence of a defect, the light pulse 
spreads in a way that strongly depends on the nonlinear 
coefficient. We arc not interested in this free propagation, 
which was studied in detail in the past fl2j ]. Ideally, we 
would like to analyze the propagation of a packet whose 
shape is, in the absence of the obstacle, independent of 
time, in order to clearly single out the influence of the 
defect on its propagation. There are different ways to 
realize this. One possibility, which is easily accessible 
experimentally, is to control the on-site energies e k by 
modulating the width of each waveguide to build a har- 



monic confining potential, e k = eo 
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defines the center of the lattice. The advantages of such 
a set up are multiple. One can shine on the lattice a light 
packet whose center is located at an arbitrary distance 
d from the bottom of the potential. As it propagates 
in the longitudinal direction, the packet will oscillate in 
the transverse direction with frequency u) . One can 
show that for 7 > the frequency uj coincides with that 
of 7 = [17[. Moreover, the shape of the packet does not 



is the reference on-site energy, and uj measures the fre- 
quency in units of normalized l/z (IBj . The site k = 



vary in time if initially it is given by the (translation) 
of the stationary ground state solution of Eq.([T]). Thus, 
with a positive nonlinearity, such a light packet oscillates 
with a frozen shape, and its velocity at the bottom of 
the potential is v = uj ■ d. This allows to control the 
transverse speed of light. 

We now include the defect at the center of the har- 
monic potential. A simple way to experimentally im- 
plement it is by a local variation of the on-site energy, 
(-k = £o + \u 2 k 2 + UoSkfi (where Uq represents the de- 
fect strength). The purpose now is to study, for different 
relative velocities v, the oscillations of the light packet 
in the presence of the defect. If the light is scattered by 
the defect, dissipative processes are induced that trans- 
form the coherent collective oscillation into disordered 
fluctuations of the light intensity. As a consequence a 
damping of the collective character of the oscillations is 
expected. On the contrary, if SF motion occurs, the light 
pulse is able to move through the defect without losing 
collectivity (e.g., without changing its global shape). It 
only creates a local intensity depletion around the defect 
[1 81 ] . By analogy with the dispersion relation of the corre- 
sponding (continuous) Gross-Pitaevsky equation, Eq.([T]) 
predicts a SF motion for velocities below a critical thresh- 
old v c which, for weak perturbations, is of the order of 
yj C-7 1^4| 2 , where \A\ 2 is the light intensity at the center 
of the light pulse. For typical experimental set ups fl9j . 
this velocity is v c w 2 • 10~ 2 , which corresponds, in the 
original units, to the ratio of the transverse speed to the 
speed of light in the photonic structure. 

Figure [5] shows, for different initial positions d, the 
first oscillations of the light packet. In the absence of an 
obstacle (Figj2^,), the packet oscillates freely with con- 
stant shape and amplitude. In presence of the defect and 
for small amplitudes (Figj^b), no damping or dissipative 
process is observed (see right column of the Figure) . The 
only manifestation of the presence of the defect on the 
light pulse is a local intensity depletion at the position of 
the defect, which is clearly visible in FigJ^b as a horizon- 
tal red line. The motion is qualitatively similar to the free 
oscillation of the light pulse, shown in Fig(5^, aside from 
a slight modification of the frequency, that can be ex- 
plained theoretically [2(J. Increasing the amplitude and 
thus the relative velocity with respect to the defect, there 
is a critical speed above which the shape of the light pulse 
is qualitatively modified as it propagates (Fig J2fc). What 
is observed is, in particular, the emission of grey soliton- 
like perturbations (dark blue trajectory), which detach 
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FIG. 2: (color online) Time evolution of a light pulse propa- 
gating through an array of waveguides modulated by a har- 
monic on-site energy, with parameters C = 100, w = 0.1, 
7 = 80, N = 300. Left column: intensity, normalized to 
the maximal one. Right column: amplitude of the oscillatory 
motion of the light packet barycenter (for the first 50 oscil- 
lations; the amplitude is normalized to the initial value). Top 
panel (a): free motion. Following panels (b)-(d): oscillation 
in presence of an additional obstacle located at the center of 
the lattice, (b): SF motion, (c): speed slightly above the 
critical velocity, (d): speed well above the critical velocity. 



from the defect and travel across the light packet with a 
non trivial dynamics. This dissipative process produces 
a damping of the oscillations. As the velocity is further 
increased (Fig(3H), the shape of the light pulse is sub- 
jected to massive deformations through phonon-like and 
solitonic emissions. The complex dynamics of the excita- 
tions signals the onset of a strong dissipative process that 
destroys the collectivity of the oscillations, thus leading 
to a strong damping. 

A quantitative way to characterize the dissipative pro- 
cess observed in Fig[2]is to numerically evaluate the flu- 
idity factor, defined as the ratio of the amplitude of the 
oscillation around some final time, to the initial ampli- 
tude, Y, k k\A k \ 2 /Y, k k\Al\ 2 , where the A^'s are the inci- 
dent amplitudes of the field. This factor varies from for 



a totally damped motion to 1 for an undamped one. We 
show in Figj3] the computation of the fluidity factor for 
different velocities and different attractive (Uq < 0) and 
repulsive (Uo > 0) defect strengths. The final time cor- 
responds to 50 free oscillation periods. The velocities are 
normalized to the perturbative critical velocity, defined 
as c s — \/2C(i, where fj, is the chemical potential of the 
incident light packet, [i = T, k A\* " {-C (A° k+1 + + 
-f\A° k \ 2 A° k + ±(wfc) 2 A°}. At low velocities, the fluidity 
factor is equal to one and the light pulse presents a per- 
fect transmission through the scattering potential. A lo- 
cal peak (resp. dip) is observed on the light intensity 
when the defect is attractive (resp. repulsive), but this 
does not affect qualitatively the dynamics of the oscil- 
lations. This demonstrates that the transverse motion 
of the light is superfluid for a well defined parameter 
range in transverse speed of light and defect strength. 
As the velocity increases a sharp transition towards a 
phase of damped dynamics is observed. This border de- 
fines the critical velocity v c . Above v c , non local dissi- 
pative excitations are allowed, and superfluidity breaks 
down. As shown in Fig(3l v c coincides with the perturba- 
tive one, c s , for weak defect strengths. However, as the 
strength increases, v c deviates from the perturbative esti- 
mate. A non-perturbative analysis is required to describe 
the threshold H3| ■ That analysis also allows to ex- 
plained the dissymmetry observed in Figj3] between pos- 
itive and negative values of Uq. The critical velocity as a 
function of the strength of the defect, shown in FigEl is 
computed by solving the equation Uo/fi = K(v c (Uq)/c s ) 
: (-8x 4 - 20x 2 + ! + (! + 8.x 2 ) 3 / 2 ) 1/2 - 



where K{x) 

The method used here to test superfluidity of light has 
a close counterpart in the physics of ultracold atoms. 
The influence of a local potential on the flow of a con- 
densate 13, l2l| or on the damping of dipole oscillations 
20l 22 , 23| became in recent years an important and ac- 



cepted experimental and theoretical tool to analyze the 
dynamics of ultracold Bose-Einstein condensates, in par- 
ticular as a test of superfluidity. 

We have considered here nonlinear optics at a purely 
classical level. In the paraxial approximation, its formal 
description is identical to the Gross-Pitaevsky equation, 
which describes a mean-field dilute condensate of atoms. 
Many open issues, related to the underlying microscopic 
quantum theory of light and to its connections with the 
phenomenon of Bose-Einstein condensation, deserve fur- 
ther investigation (l3 |. In a wider context, the physics 
described here is similar to that observed in the wave re- 
sistance of a moving disturbance at the surface of a liquid 
24j , or to the Cherenkov radiation of a charged particle 



moving through a dielectric medium [25j. 

To conclude, we have shown that, in a typical exper- 
iment, for transverse speeds of the order of 10~ 2 the 
speed of light in a self-defocusing nonlinear medium the 
light motion becomes superfluid. In contrast, superfluid- 
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FIG. 3: (color online) Fluidity factor computed after a time 
Zf = 50 -2n/uj, as a function of the normalized strength of the 
defect and transverse speed of the light pulse. Yellow regions 
correspond to SF motion, whereas damping is observed in 
the darker regions. The dashed line represents the analytical 
non-perturbative prediction of the critical velocity (see text). 
The full dots (a) — (d) indicate the parameters used in Fig[2] 
{(a): U = 0, v/c s = 0.23; (b), (c) and (d) : U /n = 0.48 and 
v/c s — 0.23, 0.76 and 1.05, respectively). 



ity does not occur in a focusing (7 < 0) medium. The 
effect described is not inherent to discrete lattice struc- 
tures, and is expected to occur in continuous media as 
well, provided the refraction index can be carefully de- 
signed. The main interest of the set up based on an 
array of waveguides is the ability to easily control the 
different parameters. Furthermore, in our effective one- 
dimensional geometry we have also identified the emis- 
sion of solitonic and phonon-like excitations as the main 
mechanisms that contribute to the breakdown of the SF 
motion above the critical velocity. In two dimensions su- 
perfluidity is also expected to occur, with a transition 
to a dissipative flow related to the emission of optical 
vortices. We believe the SF motion described here is a 
general property, which may be observed for an arbitrary 
scattering potential (not limited to a localized defect). 
The propagation in the presence of, e.g., random fluctu- 
ations of the on-site energies or of the refraction index 
is of particular interest, since randomness is inherent to 
any fabrication process. In analogy with similar recent 
studies in the physics of ultracold atoms 23[ 2(| , one may 
expect the existence of SF motion of light in presence of 
disorder. 
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